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Abstract 

We present details of our experimental and theoretical study of the components of 
the anisotropic magnetoresistance (AMR) in (Ga,Mn)As. We develop experimen- 
tal methods to yield directly the non-crystalline and crystalline AMR components 
which are then independently analyzed. These methods are used to explore the un- 
usual phenomenology of the AMR in ultra thin (5nm) (Ga,Mn)As layers and to 
demonstrate how the components of the AMR can be engineered through lithogra- 
phy induced local lattice relaxations. We expand on our previous [Phys. Rev. Lett. 
99, 147207 (2007)] theoretical analysis and numerical calculations to present a sim- 
plified analytical model for the origin of the non-crystalline AMR. We find that the 
sign of the non-crystalline AMR is determined by the form of spin-orbit coupling 
in the host band and by the relative strengths of the non-magnetic and magnetic 
contributions to the impurity potential. 
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1 Introduction 



The Anisotropic Magnetoresistance (AMR) in ferromagnetic metals has been 
the subject of many studies since it was discovered more than a century ago[lJ. 
The effect describes the change of the electrical resistance in response to a 
change in the orientation of the magnetization. Phenomenologically, AMR has 
a non-crystalline component, arising from the lower symmetry for a specific 
current direction, and crystalline components arising from the crystal symme- 
tries [2|3] . In ferromagnetic metals the AMR components have been extracted 
indirectly from experimental data by fitting the total angular dependences 
[3]. The coefficients can be obtained by numerical ab initio transport calcula- 
tions [4j, but these have no clear connection to the standard physical model of 
transport arising from spin dependent scattering of current carrying low mass 
s-states into heavy- mass d-states [5]. 

In the dilute magnetic semiconductor (Ga,Mn)As, it has been observed that 
the AMR consists of a non-crystalline component of the opposite sign (com- 
pared to most ferromagnetic metals) and typically much weaker crystalline 
terms reflecting the underlying magneto crystalline anisotropies piTfSI^ITIIfTTjlT! 
In our recent work [13] we developed experimental methods to yield directly 
the non-crystalline and crystalline components of the AMR in (Ga,Mn)As. 
Numerical microscopic calculations were employed to calculate the AMR com- 
ponents, achieving good agreement with experimental results. An analytical 
model was introduced to explain the physical origin of the AMR terms, in 
particular to show that the sign of the non-crystalline term is determined by 
the ratio of the magnetic and non-magnetic components of the scattering po- 
tential. Here we expand our previous work to include details of that analytical 
model and a detailed explanation of how the components of the AMR can be 
separately extracted from the experimental data. 



2 Experiment 

25nm and 5nm Gao.95Mno.05As films were grown by low temperature molecular 
beam epitaxy onto a 330nm HT-GaAs buffer layer grown on a semi-insulating 
GaAs(OOl) substrate. Standard photolithography techniques were used to fab- 
ricate two sets of Hall bars of width 45yum with voltage probes separated by 
285/im, with the current along the [100], [010], [110] and [110] directions. The 
longitudinal resistance R^^^^ and the Hall resistance R^j; were measured using 
four probe techniques with a DC current of 1/iA. We also measured devices 
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Fig. 1. (a) Measured (at 4.2 K) longitudinal and (b) transverse AMR for as-grown 
25nm thick Gao.95Mno.05As film as a function of the angle between magnetization 
in the plane of the film and the current. The legend shows the direction of the 
current. The y-axes show lS.pl pav shifted so that the minimum is at zero, to show 
the symmetries present in the data more clearly, (c) Definition of the angles referred 
to in the text, (d) Corbino geometry. 

with a Corbino geometry (Fig. 1(d)) in which the current flows radially in 
the plane of the material. All measurements were made on as-grown material 
(i.e. without annealing). Magnetoresistances were measured with a saturating 
magnetic field of IT applied in the plane of the device, i.e., in the pure AMR 
geometry with zero (antisymmetric) Hall signal. Magnetometry measurements 
were carried out using a commercial Quantum Design SQUID magnetometer. 

Figures 1(a) and (b) show the fractional change in the resistivity, ^p^xl Pav 
and ^pxyl pav i^Pxx = Pxx — pav) for the Hall bars fabricated from the 25nm 
film for each current direction. Here pav is the average value of the longitudi- 
nal resistance as the magnetic field is rotated through 360°. In the following 
analysis we will show that this data can be decomposed into the form given 
by equations ([1]) and ([2]). In the subsequent sections we aim to identify the 
origin of the particular contributions. 

The terms allowed by symmetry are obtained by extending the standard phe- 
nomenology [2], to systems with cubic [100] plus uniaxial [110] anisotropy: 
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Ci COS 2(j) + Cu cos 2ij + Cc cos i^j + Ci^c cos(4^/^ - 20) . (1) 
Ci sin 20 - Ci^c sin(4^ - 20) . (2) 

is the angle between the magnetization unit vector M and the current J, and 
■0 the angle between M and the [110] crystal direction. The four contributions 
are the non-crystalline term, the lowest order uniaxial and cubic crystalline 
terms, and a crossed non- crystalline/ crystalline term. We have omitted the 
higher order terms in these expressions as these are found to be negligibly 
small in our devices. The purely crystalline terms are excluded by symmetry 
for the transverse AMR and it is clear that lS.pxyl Pav and lS.pxxl Pav are not 
independent. The well known formulae for isotropic materials (see e.g. [8j) 
correspond to C/ 7^ and all other coefficients set to zero in equations ffTP]) . 
Such formulae result also from equations (fT][2!) averaged over -0 with held 
constant, which corresponds to a polycrystal with randomly oriented grains. 

It is not possible to obtain each term separately by fitting equations ([1]) and 
([2]) to the full (^z^) angular data from one Hall bar. In a [110] oriented Hall bar, 
for instance, = -0 + |, the equations reduce to /S.pxx/ Pav = {.~Ci + Cu — 
C/,c) cos 2?/) + Cc cos Alp and Pxy/ Pav = —Ci + C/^c* sin 2-0. Ci and Ci^c could 
be determined only if Cjj = 0. However, using different Hall bar orientations, 
the anisotropy constants appear in other combinations so that they can all 
be extracted individually. In the present case, two Hall bar orientations (e.g. 
[110] and [110]) suffice. Measurements with other orientations are used as 
consistency checks. 

Figure 2(a) shows some examples of how the individual AMR components can 
be extracted: the longitudinal resistivities in the [010] and [100] directions are 
subtracted to give a simple sin 2-0 signal (cf. Eq. [T]) with amplitude 2{Cj—Cj^c) 
and this is the same as the difference of transverse resistivities in the [110] and 
[110] samples. A cos 2-0 signal with amplitude 2{Ci + Ci^c) can be obtained 
through other combinations of the data. The close agreement obtained by 
using different combinations to extract the same coefficients highlights the 
consistency between the data and our phenomenological analysis, and also 
attests to the homogeneity and quality of the MBE grown GaMnAs wafer. 

The purely crystalline terms can be extracted from the combination of the 
measurements from two orthogonal bars. This is shown in Fig. 2(b) for our 
data. In principle, the same signals could be extracted directly from a single 
experimental measurement by appropriate wiring of an L-shaped Hall bar. 
Alternatively, the same coefficients can be extracted directly from the data 
obtained from the Corbino geometry. Fig. 2(b). As the current flows in all 
radial directions, the non- crystalline terms are averaged to zero[Tl]. It would 
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Fig. 2. Combinations of the data from Fig.l used to extract the AMR components. 
The notation xx[100] is used to denote Ap^x/pav for current along the [100] direction 
etc. 



be problematic to perform such measurements with ferromagnetic metals be- 
cause, in order to define the Corbino geometry it is required that the source- 
drain resistance be large compared to the contact resistances. This condition 
is met with (Ga,Mn)As. Again, we obtain excellent agreement for the same 
coefficients extracted from different combinations of devices. 



Figure 3(a) shows the coefficients Cj^c, and extracted from the Hall 
bar and Corbino disk data over the whole range up to the Curie temperature 
(80K). Note that the uniaxial crystalline term, Cjy, becomes the dominant 
term for T^30K. This correlates with the uniaxial component of the magnetic 
anisotropy which dominates for T^SOK as observed by SQUID magnetometry 
measurements (Figure 3(b)). Our work shows that in (Ga,Mn)As ferromag- 
nets, the symmetry breaking mechanism behind this previously reported [13] 
uniaxial magneto-crystalline anisotropy in the magnetization also contributes 
to the AMR. 



In our previous work ^13j we used the techniques explained above to extract 
the individual AMR components for ultra thin (5nm) Gao.95Mno.05As films for 
which the AMR (Fig. 4(a)) is very different from that observed in the 25nm 
film. For the 5nm film, the uniaxial crystalline term dominates the AMR as 
shown in the Fig. 4(b) for the device in the Corbino geometry. These techniques 
were also used to demonstrate how the crystalline terms can be tuned by 
the use of lithographic patterning to induce an additional uniaxial anisotropy 
in very narrow Hall bars. It has been shown P^fTU] that the fabrication of 
narrow bars allows the in-plane compressive strain in the (Ga,Mn)As film to 
relax in the direction along the width of the Hall bar and this can lead to 
an additional uniaxial component in the magnetocrystalline anisotropy for 
bars with widths on the order of l/xm or smaller. Figs. 5(a) and (b) show 
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Fig. 3. (a) Temperature dependence of the crystalline terms extracted from the Hall 
bars and Corbino devices, (b) Remnant moment of 25nm Gao.95Mno.05As film mea- 
sured by SQUID magnetometer along different crystalline directions after cooling 
in a field of O.IT. 




the AMR of 45/um wide bars and 1/xm wide bars fabricated from nominally 
identical 25nm Gao.95Mno.05As wafers. For the 45/im bars, the cubic crystalline 
symmetry leads to the AMR along [100] and [010] being larger than along [110] 
and [110]. For the narrow bars we observe the opposite relationship, consistent 
with the addition of an extra uniaxial component which adds 0.8% to the AMR 
when current is along [110] and [110] and subtracts 0.4% when the current is 
along [100] and [010]. 



It is anticipated that the experimental techniques developed here will be useful 
in studying the magnetotransport coefficients in other magnetic materials and 
nanostructures. 
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Fig. 5. (a) AMR for macroscopic Hall bars and (b) narrow (1/im wide) Hall bars 
fabricated from 25nm Gao.95Mno.05As films. T=4.2K. 

3 Prom the full Boltzmann theory simulations to a simple analyt- 
ical model of the AMR 

A six-band k ■ p description of the GaAs host valence band combined with the 
kinetic-exchange model of the coupling to the local Mnca ci^-moments [20] can 
provide input for the Boltzmann equation producing the conductivity tensor 
a. In our previous work [13], we expanded the studies of Jungwirth et al. 
[23f7] by performing numerical calculations of the axx for arbitrary in-plane 
directions of M (not just M\\I and M ± /). The results of the calculations 
(see Fig. [6]) are in semi-quantitative agreement with the experimental AMRs 
for the 25nm Gao.95Mno.05 As films shown in Fig.[T](a,b). The zero-temperature 
model of a 15% compensated material [24\ shows that both the Ap^x/Pav and 
Apxy/ Pav are dominated by the non- crystalline cos 20 and sin 20 part, being of 
the order of several percent and negative [Cj = —2.2%). The crystalline terms 
are an order of magnitude smaller {Cu = —0.2%, \Cc\ < 0.1%, Cj^c = 0.4%). 
These calculations assumed a growth strain cq = —0.3% and a uniaxial strain 
[21] = -0.01%. 

Previously, we noted that the crystalline terms arise from the warping of 
the valence band [13]. In the rest of this paper we concentrate only on the 
non- crystalline component of the AMR (Ci) which in our model originates 
from anisotropic scattering of spin-orbit coupled holes on MuQa impurities 
containing polarized local moments. 

The six-band Kohn-Luttinger Hamiltonian Hkl is parameterized by the Lut- 
tinger parameters 71, 72, 73 and SO splitting Aso while the exchange splitting 
strength is characterized by the constant Jpd [20] . The total Hamiltonian then 
reads 

H = Hkl + Jpd E Si ■ s^6{f, - Rj) , (3) 
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Fig. 6. Calculated angular dependence of the AMR (see the text for system param- 
eters) . 

where Sj and Si denote the Mn and hole spins placed at rj and Rj. The 
mean-field treatment of the exchange splitting (second term in Eq. ([3])) leads 
to a term formally equivalent to an effective Zeeman splitting of strength 
h = JpdNunSun, where A^Mn is the Mn concentration and ^Mn = | for the five 
Mn d-electrons. In the following we consider the spherical approximation to 
Hkl (72 = 73) and the four-band model (A50 ~^ 00). The total Hamiltonian 
can be then rewritten as 

H ={-fi + §72)^2^72771 - 272(fc ■ Jfh^/2m + hM ■ s. (4) 

Explicit expressions of the matrices for the total angular momentum J and 
hole spin s can be found in Refs. [26] and [22]. With this simplified band 
structure we will now estimate the AMR assuming that the scattering takes 
place only on Mn ions substituting Ga which is the dominant mechanism. 

Numerically, the light hole bands are found to give only small contribution 
to the total conductivity as shown in Fig. ([7]), hence we consider only the 
heavy hole bands. This is a substantial simplification since the Hamiltonian 
(jl]) restricted to the heavy-hole subspace j'^ = ±| can be diagonalized analyt- 
ically (see Appendix A). By denoting the angle between k and M as 0^ the 
eigenenergies read 

^1,2 = (7i- 272)^'A;V2m±/7cos0^. (5) 

In the presence of the exchange field h the two originally identical heavy-hole 
Fermi spheres transform into (typically weakly) distorted '+' and '— ' spheres 
displaced by A/c = :k{2m^/hY^'^h/Ay/^F along M. The eigenstates \hh + k), 
\hh — k) are independent of h and correspond to perfectly isotropic radial spin 
textures: expectation value of s in the state \hh + k) {\hh — k)) is a vector 
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Fig. 7. Contribution of particular bands to the calculated total bulk conductivity in 
uncompensated GaAs. (a) Absolute values, (b) relative to the total conductivity. 

Magnetisation 
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Fig. 8. Heavy hole bands. Spin texture without and with exchange splitting. 

of length i parallel to k and pointing outwards (inwards), Fig. [HI We remark 
that Eq. ([5]) and \hh±k) remain unchanged even when the Hamiltonian ([5]) 
is extended to include a finite Ago (see Appendix A) . 

In the lowest order Born approximation, the transport relaxation time is given 

by 



(Pk' 27r 



r{k) J (27r)3 h 



NuJiEik) - Eik'))\Mr,,\\l -k-k'). 



(6) 



The last factor corresponds to vertex corrections in the Kubo formalism and 
k- k' is the cosine of the angle between k and k'. The scattering matrix element 
between initial \zj^,) and final l^;^) states is 

Mg, = {z^,\V{k-k')\z^,,) = V{\k-k'\){z^,\z^,,) (7) 
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for the isotropic (e.g. screened Coulomb) scattering off a charged ion, and 

M|:, = {z^\{h/NMn)M ■ s\z^^,) = Jj,,SMn{z^\M ■ s\z,^,) (8) 

for the scattering off a magnetic moment. These two add up coherently for 
substitutional Mn so that Mrp = + Mg,. 

KK kk' kk' 

To simplify the further qualitative analytical discussion we replace the long- 
range Coulomb potential with an effective 5-function potential and denote the 
ratio of this non- magnetic potential and of JpdSMn as a. Our relaxation rates 
will be evaluated from the following scattering operator 

^kk' + ■ (9) 



We discarded an overall prefactor, since when only one source of scattering 
is present, this prefactor will cancel in the expressions for relative changes of 
conductivity needed to obtain the AMR. 

Equation ^ will now be used to calculate the scattering rates for two special 
values of fc, details of that calculation are given in Appendix B. The first one 
is k\\M = e^. The corresponding initial states for the scattering, \hh±k) will 
be abbreviated as |0±), see Fig. [HI Integrands in Eq. can be simplified 
considerably noting that s^^lOi) = ±i|0±). We obtain 



_= / d^k' l\{hh + k'\0±)Wa±-y+ (10) 
rf J [ 2 

FS 

+ \{hh - k'\0±)\^{a ± ^)2| {1-k-k'). 

The integration variable is taken as dimensionless and the result should be 
multiplied by a factor R = 2mkF/h? ■ 2'K/h ■ Nmu ■ {JpdSMn}'^ to get the real 
inverse scattering times based on Eqs. (ITfHjl with V{k — k') = aJpaSun- By 
integrating Eq. ffTOj) we obtain. 



1/r^ = 27r(a - \ f 1/rf = 2n{a + |)2 . 
In contrast, the other pair of initial states |1±) with ± M (see Fig. [8]), gives 



L = J d^k' {\{hh + k'\l±)\^a^ + \{hh + k'\s,\l±)\^+ 

'^^ FS 

\{hh - ^'|1±)| V + \{hh~ k'\s^\l±)\^ + } {1-k-k'). (11) 
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The integrals not involving Sx give again 27r by the virtue of symmetry (cf. 
Eq. ( ITOl) and Fig. [8]) while the remaining two terms have to be evaluated 
independently to give 



l/rl = Ijrl = 2vr(a^ + 1) (12) 

Conductivities for M\\I and M I can be estimated based on these relaxation 
rates as a\\ = P{t^ + t^) and a± = P{t^ + r^), i.e. 



where P = e^Vp/S ■ gR/27r, vp is the Fermi velocity and g is the density of 
states at the Fermi wavevector kp- For heavy holes, vp = (71 — 2'~f2)^kF/m 
and g = 2mkF/h'^{ii - 272). Eq. ^ leads to 



ail _ {a^ + ^){a^ + \) 
(a2 - \f 



(14) 



presented in Ref. [13]. The AMR, using the notation of the experimental sec- 
tion then reads 

AMR = . P--(^\\I) - ^ I) ^ , a||-a± ^ 20a2 - 1 



When the magnetic term in the impurity potential is much larger than the 
non-magnetic term (a <^ 1) one expects a\\ < a± (positive AMR, as is usually 
observed in metallic ferromagnets). However, the sign of the non-crystalline 
AMR reverses at a relatively weak non-magnetic potential {a = 1/a/20 in the 
model), its magnitude is then maximized when the two terms are comparable 
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{a = 1/2), and, for this mechanism, it vanishes when the magnetic term is 
much weaker than the non-magnetic term {a oo). Note that the large 
magnitude of the AMR for a = 1/2 due to l/cry = is an artefact of the 
special form of the simplified scattering operator (9) where the magnetic and 
non- magnetic part depend in exactly the same way on k (they are constant, 
i.e. they both correspond to point-like scatterers). 

Physically, carriers moving along M, i.e. with s* parallel or antiparallel to M, 
experience the strongest scattering potential among all Fermi surface states 
when a = 0, giving cry < a±. When the non-magnetic potential is present, 
however, it can more efficiently cancel the magnetic term for carriers moving 
along M, and for relatively small a the sign of AMR flips. Since a < 1/ v^20 is 
unrealistic for the magnetic acceptor Mn in GaAs [23|fT7] we obtain ay > a±, 
consistent with experiment. 



4 Summary 

A phenomenological framework for the AMR in the zinc-blende crystalline 
environment of GaMnAs was used to analyse experimental data from (i) bulk 
material, (ii) thin layers, and (iii) samples with lithographically manipulated 
strain. While leaving the crystalline components aside for a further theoretical 
study, a qualitative analytical model was presented for the non-crystalline 
AMR. The model is based on the anisotropic relaxation times in the heavy hole 
bands due to a combined magnetic/non-magnetic scattering on Mn impurities. 
It offers an explanation of the sign of the non-crystalline AMR in GaMnAs 
which is opposite to most of the conventional metal ferromagnets. 

We acknowledge support from EU Grant IST-015728, from UK Grant GR/S81407/01, 
from CR Grants 202/05/0575, 202/04/1519, FON/06/E002, AV0Z1010052, 
KAN400100652, and LC510, from ONR Grant N000140610122, and from SWAN. 
J. Sinova is a Cottrell Scholar of Research Corporation. 



A Some basic algebra with Kohn— Luttinger Hamiltonian 

Equation and explicit form of the heavy-hole eigenvectors \hh±k) will be 
derived in this Appendix. As we will refer to the spherical approximation with 
Mlle^:, the main body of calculations will be carried out in the kx,ky plane 
keeping in mind the rotational symmetry of the problem around ix. For the 
sake of completeness, results for general k^ will be given at the end of 
each part without derivation. 
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The starting point is Hkl as in Eq. (A8) of Abolfath et aZ.|22 
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Going into the kx,ky plane by putting k = /c(cos</), sin0, 0), the off-diagonal 
elements become 



6 = 0, c = [72 cos 20 — 'j73 sin 20] , d 

2m 

and the diagonal elements are independent of 



^/2h^k^ 



2m 



-72 



2m 



■(71 + 72), i^zfe = ^^(71 - 72) 



2i,2 



if. 



2m 



■71 + A 



so 



It is convenient to factor out h^k'^ / (2m) and to introduce 5 = A5o/(72^^/2m) ■ 
k~\ 

The spherical approximation relies in setting 72 = 73 and allows for the sim- 
plification 

c/{h^k'^/{2m)) = V372exp(-2i0) . 
In this approximation and even with general k^ 7^ 0, i.e. 



k = k{sm 6 cos 0, sin 6 sin 0, cos ( 



(A.2) 



the spectrum of the Hamiltonian (lA.ip does not depend on 6, 0. Below, this 
statement is demonstrated explicitly for 9 = 7i/2. 



A . 1 Spectrum 



When kz = 0, rows/ columns 3,4,5 of H^l are totally decoupled from rows/ columns 
1,2,6. The 6x6 problem is reduced to two independent and equivalent 3x3 
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Fig. A.l. Dispersions of the heavy holes, hght holes and the split-off band in 
the spherical approximation. Horizontal dashed line indicates the Fermi energy of 
Ep =100meV. 



problems. The latter (lines 1,2,6) is 

/ 

A+ = H/{h^k''/2m) = 



7i + 72 -c c\/2 
-c* 7i - 72 -V272 
V2c* -v^72 7i + I2S 



(A.3) 



Matrix depends on only via c. In the spherical approximation, cc* = 
|cp = 37I is ^-independent. The eigenvalues of are calculated from det(y4+ — 



A) = 0. This determinant will depend on c only via the combination cc* 
hence it is 0-independent for 72 = 73- 



and 



Explicitly, Fig. IA.lt the eigenvalues of A+, multiplied back by h'^k'^/2m, are 



Ei = ^ [71 - 272] 



E2 
E3 ■■ 



2m 



2m 
2m 



li + l2--{S + 2- VS^ - 45 + 36) 



(A.4) 



71 + 72 ■ 2 (-^ + 2 + -45-1- 36) 



They correspond to heavy holes, light holes and the split-off band, the last two 
depend on k via 6, but none of them depends on (p and hence give circular FS 
sections. The light-hole energy is [71 + 272]^^fc^/(2m) in the limit Aso 



00. 



Analysis of the complete Hamiltonian Hkl is completed by constructing ma- 
trix A_ from lines 3,4,5. It is not identical to A^, but it has an identical 
spectrum to A+. Thus, the Fermi surfaces of H^l in the k^ = sections 
comprise of, largest to smallest, two coincident circles for heavy holes with 
Ep = (71 — 2'j2)fT''^kp/2m, which is independent of Aso, two coincident circles 
for light holes and two coincident circles for split-off bands. 
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A. 2 Eigenvectors 



Even though the spectrum of the Hamiltonians (lA.ll) or (]A.3P does not depend 
on 0, 9, the eigenvectors do. 

The lowest-energy {Ei = 7i — 272) eigenvector of is independent of 6. 
Together with its counterpart from A_ they represent the degenerate heavy- 
hole states of Hkl regardless of Aso- 

vi+ = i(e-2^^y3,0,0,0,0f , vi- = ^(0, 0,73,6^^^0, of . 



The component notation refers to the basis (A5) of Ref . p2] , the superscript 
denotes the transposition. For a general k given by Eq. (1A.2P the heavy-hole 

states are 



I sin 9, VS sin 9, 0, -e^"^ cos 9, 0,0)^, 



vi^ = ^ {e"'"'' cos9,0,V3sm9, e"^'^ sin ^, 0, 0)^ . 



A. 3 Perturbation theory with kinetic exchange 



The p-d kinetic exchange has the form of an effective Zeeman field in the 
magnetisation direction which we choose to be x, Fig. [HI The corresponding 
operator (with respect to the first four basis vectors of (A5) in Ref. [22]), that 
is the last term in Eq. (jlj), is 



pd 



6 



/ \ 

2 73 

73 2 

y 73 ] 



(A.5) 



The matrix of Hp^/h in the basis of ^1+, Vi- and its diagonalization results 
are the following 



Hpd/h 



1 + e"^'^ 
1 + e-2^<^ 



f cos (/) , 



E_ = — i cos 



= 2-1/2(1, -e-**)^ 



In the limit of small h ('degenerate- level perturbation calculus') the degenerate 
heavy-hole bands of Hkl become split in energy by h cos 0. Recast into the k- 
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plane, the two coincident Fermi circles become displaced. Their wavefunctions 
are 



\hh+, k) = ^(e^^*^, VS, VSe-'''', e'^, 0, 0)^ , (A.6) 
\hh-, k) = ^(e-'*^ v^, -v^e-'^ -e^^ 0, 0)^ . 



For general k^ ^ Q flA.2p the energies are 



E± = ±- cos (j) sin 9 , 
as stated in Eq. ([5]), and the eigenvectors 

\hh+, k) = {e-^'^C\ VSCS^ VSe-'^C^S, e'^S\ 0, 0)^ , (A.7) 
\hh-, k) = (e-^^'^S^ V3SC\ -VSe-'^^S^C, -e^'^C^ 0, 0)^ 

with C = cosd/2 and 5* = sin6'/2. 



B Integrals of overlaps, Eqs. ( llOlllip 

The overlaps of eigenvectors in Eq. (1A.6P are 

{hh + k\hh + k') = -[{e^'(^-^') + 3) + e*('/'-</'')(e-2i(</'-0') + 3)] ^ 
8 

{hh -k\hh + k') = -[{e^'^'t'-*') + 3) - e^(</'-*')(e-2K</'-<^') + 3)] . 
8 

Note that these overlaps for J = \ spinors are not the same as for J = \ 
spinors |a ± A;) = (Ijie*"^)"^ even though the spin textures of \hh ± k) and 
|a ± fc) are the same, Fig. [HI 



The two summands of integral in Eq. (iTOll . 

(fk'\{hh + k'\Q±)Wl-k-k') (B.l) 



FS 



will be integrated with spherical coordinates 0', $' with 'north pole' 0' = in 
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fe'llca;. The polar angle ($') integration gives a factor of 27r and we obtain 

— / d(j)' sin (J)' ■ 2 [10 ± 15 cos 0' + 6 cos 20' ± cos 30'] (1 - cos 0') 
64 J 



For the + sign we get 47r(6 + |)/64 and for — it is 47r(26 — |)/64 giving the 
total of 2tc. The integrals with {hh — k'\ in Eq. fIB.ll) are analogous. 



Concerning the integral in Eq. (fTTi) . we first have to show that it is indeed 
equal to the scattering rate 



^ = J d^k' [\{hh + k'\a + s^\l±)\^+ 
\{hh + k'\a + s^\l±)\'^} {1 - k ■ k') . 



FS 



In order to show this, it suffices to demonstrate that 



2lT 



J d<^>'{hh + k'\a\l+){l + \s^\hh + k') = 



which holds by the virtue of {hh + A;'|s^|l+) oc e**' and {hh + k'\l+) oc e^**' 



Proceeding with Eq. ffTTl) . we will use 

{hh + k\s^\hh + k') = he'^ + e-^'^')[l + cos(0 - 0')] 

8 

{hh — k\sx\hh — k') = —{hh + k\sx\hh + k') 

{hh - k\s.x\hh + k') = --{e'^ - e"^')[l - cos(0 - 0')] 

8 

{hh + k\sx\hh — k') = —{hh — k\sx\hh + k') 

Here, the proper choice of k {\\ey) amounts to putting = | and the integral 
for scattering from |H-) to the + band is 

tt/2 

— [ rf0'cos0'[l + sin0'l22(l-sin0') ■ (1 -sin0') = — • — 
64 i ^ ^ ^' ^ 32 120 

-7r/2 

while the other integral {—>■—) 
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7r/2 

^ / #'cos0'[l - sin 0f 2(1 + sin 00 ■ (1 - sin0') = ^ " " 3^) ■ 

-7r/2 

Their sum is 27r^ which completes the proof of Eq. (fT2ll . 
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